Constraints on a Stochastic Background of Primordial Magnetic Fields with WMAP 

and South Pole Telescope data 
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We constrain a stochastic background of primordial magnetic fields (PMF) by its contribution to 
the cosmic microwave background (CMB) anisotropy angular power spectrum with the combination 
of WMAP 7 year and South Pole Telescope (SPT) data. The contamination in the SPT data by 
unresolved point sources and by the Sunyaev Zeldovich (SZ) effect due to galaxy clusters has been 
taken into account as modelled by the SPT collaboration. With this combination of WMAP 7 
yr and SPT data, we constrain the amplitude Gaussian smoothed over 1 Mpc scale of a stochatic 
background of non-helical PMF to Bim P c < 3.5 nG at 95% confidence level. Our analysis shows that 
SPT data up to £ = 3000 bring an improvement of almost a factor two with respect to results with 
previous CMB high-^ data. We then discuss the forecasted impact from unresolved points sources 
and SZ effect for Planck capabilities in constraining PMF. 



I. INTRODUCTION 

The origin of the large scale magnetic fields observed in 
galaxies and clusters of galaxies is an open issue of great 
importance in modern astrophysics. The dynamo effect 
provides a mechanism to explain the observed magnetic 
fields associated to galaxies, whereas those associated to 
clusters may be generated by the amplification of small 
fields due to gravitational compression. Both these mech- 
anisms might require an initial magnetic seed, although 
with different amplitude and different correlation length 
for galaxies and clusters. 

In addition to observations of unusually strong mag- 
netic field in galaxies at high redshift p], i], recently 
high energy observations from FERMI added a signifi- 
cant piece to the puzzle of extragalactic magnetic fields. 
The presence of diffuse extra-galactic magnetic fields in 
intercluster voids is proposed to justify the low flux of 
GeV photons from distant TeV blazars 0, 0, i, H [| 
These data pose a lower limit of the order of 10~ 18 — 10 
Gauss to the amplitude of extragalactic magnetic fields. 
Such extra-galactic magnetic fields in voids are compat- 
ible either with a primordial origin or an efficient trans- 
port mechanism for magnetic fields in galaxies. There- 
fore, the high energy observations by FERMI are not in 
contraddiction with a primordial origin of extra-galactic 
magnetic fields. 

Current CMB anisotropy measurements lead to upper 
limits on the amplitude of a stochastic background of 
primordial magnetic fields generated before nucleosyn- 
thesis Hi, Ea El El (see Ref. [lj] for a treatment 
claiming much tighter constraints). Indeed, a stochastic 
background of PMF generates all types of mag netized 
linear perturbations [pj, EH : scalar [8|, E3, EH, EH E3, EM 
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|19|, |2fJ, vector [2JJ, |22|, |23( and tensor [22|, |2J, |25| and all 
these contribute to the CMB anisotropy pattern in tem- 
perature and polarization. PMF modelled as a fully in- 
homogeneous component have also a fully non-Gaussian 
contribution to CMB anisotropies with a non zero higher 
statistical moments, which can be used as useful probes, 
such as the magnetized bispectrum (2|| [27j and the mag- 
netized trispectrum [2S| . 

In our previous works [l(3, EH Eil we have refined the 
computation of magnetized CMB anisotropies keeping 
into account only relativistic degrees of freedom and the 
correlators for the Fourier transforms of the EMT in pres- 
ence of a sharp cut-off which mimics the damping scale 
due to viscosity 0, E3 ■ In the most recent work [l(| we 
have computed the constraints coming from CMB data 
by WMAP7 in combination with data from ACBAR [29jj . 
QUaD_[3fJ and BICEP [3lQ] updating previous investiga- 
tions mi Ei- 
in this work we use the publicly available CMB 
anisotropy data at high multipoles as those from the 
South Pole Telescope (SPT) (32|, |33[ to further constrain 
a stochastic background of PMF. See Ref. [34[ for an 
analogous constraint on another important high-^? contri- 
bution to CMB anisotropies such as cosmic strings. Con- 
straints on PMF from anisotropies at high multipoles, 
£ ~ 3000, are not a straightforward extension of what 
derived at smaller angular scales: Right at those angu- 
lar scales where the Silk damping suppresses the stan- 
dard primary CMB fluctuations, the dominant contribu- 
tion of PMF and other non-standard cosmological models 
(such as cosmic strings) is polluted by extragalactic con- 
tamination and secondary anisotropies, such as Sunyaev- 
Zeldovich. In order to fully exploit small scale CMB data 
to constrain physics which produce CMB contributions 
which are not erased by the Silk damping, it is necessary 
to model carefully the residual foreground contamination 
to the angular power spectrum. 

The paper is organized as follows. In section II we sum- 
marize the general properties of PMF and the methodol- 
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ogy based on Refs. EH [H| to compute their contri- 
bution to CMB anisotropics. In Section III we briefly in- 
troduce the contamination of foreground and secondary 
residuals at high multipoles in the perspective of SPT 
data. In Section IV we present and discuss the results 
obtained with WMAP 7 + SPT data. In Section V we 
discuss the impact of foreground and secondary residuals 
at high multipoles on the Planck capabilities to con- 
strain a stochastic background of PMF and in section VI 
we draw our conclusions. 



II. STOCHASTIC BACKGROUND OF PMF 
AND MAGNETIZED CMB ANISOTROPIES 



In this section we summarize the methodology used in 
our previous papers to compute the PMF contribution to 
CMB anisotropies. We model a stochastic background of 
PMF as a fully inhomogeneous component with a power- 
law power spectrum P B {k) = A k riB , where A is the am- 
plitude and n_B is the spectral index. Our convention for 
the Fourier transform of the two point correlation func- 
tion for a stochastic background is: 



(Bi(k)fl;(k')) = (27r) 3 <5(k - k')(^ - kh)^- (1) 



where n B > —3. We assume the MHD limit in which 
B(x, r) = i?(x)/a(r) 2 with a(r) being the scale factor 
(normalized to ao = 1 today) and r the conformal time. 
As convention, we use the amplitude of the magnetic 
fields smoothed over A: 



space [Hill: 

|n(y)(fc)|2 = 51^ J dpP B (p)P B (\\L-p\) 

x [(l + /3 2 )(l- 7 2 )+ 7 /3( Ai - 7 /3)] (5) 

|n(T)(fc)|2 = 5i^i/ pPB(p)Ps(|k " Pl) 
x (l + 2 7 2 + 7 2 /3 2 ), (6) 

IMfc)| 2 = ^2^2 ^dpP B (p)P B (\k-p\) 

x (l + M 2 +4 7 /?(7/3-M)), (7) 

where fi = p • (k — p)/|k — p|, 7 = k • p, 
(3 = k • (k — p)/|k — p and SI denotes the volume 
with p < kr>. 

The analytical exact results for the PMF EMT spec- 
tra were given for specific values of n B in 0EJ. As 111 
our previous work [10j we use accurate approximations 
for the power spectra of pe , Lb , II V which are created to 
calculate the PMF contribution to CMB anisotropies in 
a continuous range of n B . We use the initial conditions 
for cosmological fluctuations as given in [l^, [l4|. For 
scalar perturbations we consider the compensated mode 
described in [l(| EH Ell- The scalar magnetized per- 
turbations are the dominant PMF contribution to CMB 
anisotropies on large and intermediate angular scales, 
whereas the vector magnetized perturbations represent 
the dominant PMF contribution on small angular scales. 
On these scales the primary CMB is suppressed by the 
Silk damping, making the vector magnetic mode the 
dominant contribution to the angular power spectrum 
as shown in Fig. [1] To constrain the PMF amplitude we 
neglect the tensor contribution, since it is subdominant 
with respect to scalar and vector ones E3 ■ 



_ [<*> dkk 2 w A f n B + 3 \ 

(2) 

PMF survive the Silk damping but are damped on 
smaller scales by radiation viscosity [U [35|. We mod- 
elled this damping with a sharp cut off in the power spec- 
trum at the scale [2l|, |35[ : 



(3) 

A stochastic background of PMF acts as a fully in- 
homogeneous source to metric scalar, vector and tensor 
perturbations. The source terms are given by the PMF 
energy momentum tensor and Lorentz Force in Fourier 



III. ASTROPHYSICAL CONTAMINATION OF 
CMB DATA ON SMALL ANGULAR SCALES 

It is well understood and proved that CMB data are 
a fundamental tool to constrain a stochastic background 
of PMF. Considering the nature of their impact on CMB 
angular power spectrum, which is mainly on small angu- 
lar scales and dominated by the vector contribution, it 
is obvious that higher the resolution of the data tighter 
should be the constraints on PMF. But data on small 
^angular scales are also affected by contamination from 
astrophysical sources. In particular for SPT data the 
astrophysical contamination is given by residual extra- 
galactic point sources and cluster of galaxies. The resid- 
ual point sources contribution strongly depends on the 
frequency considered: frequencies lower than 90-100 GHz 
are mainly dominated by radiogalaxies whereas frequen- 
cies greater than 200 GHz are dominated by infrared 
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galaxies, the frequencies inbetween are affected by both 
contributions [36J. Both radio and infrared galaxies con- 
tribute with a Poissonian term which is due to their 
random distribution in the sky. The Poissonian term 
is simply given by a flat angular power spectrum whose 
amplitude is determined by the source number counts 
integrated in flux densities and the flux density detec- 
tion threshold [37|, |38|, |39|. Infrared galaxies together 
with the Poissonian term contributes also with a clus- 
tering term. The clustering is much more complex than 
the Poissonian term, it grows with the frequency and be- 
comes the dominant component for higher frequencies. It 
depends on the properties of infrared galaxies, on their 
redshift, on the bias and on the cosmological model, and 
can be modelled in different ways with increasing com- 
plexity [H, i3, HI SI SI • The galaxy clusters contribute 
with the Sunyaev-Zeldovich effect (SZ) which can be di- 
vided into thermal [43j] and kinetic [44j contributions. In 
the case of SPT data both contributions have been con- 
sidered in a single SZ term |32| . 

For our analysis we use the templates given by the 
SPT collaboration for the 150 GHz data [Hi, S|. The 
templates are characterized by one amplitude parameter 
each, therefore we account for three new parameters in 
the analysis. 

In the third panel of Fig. Q]we show the shapes of the 
the three astrophysical contributions, Poissonian, clus- 
tering and SZ effect. The flat shape of the clustering 
term on large scales is given in the template provided 
by the SPT collaboration which considers a shape varia- 
tion only on multipoles where the clustering contribution 
is relevant, assuming a flat contribution elsewhere, con- 
sidering that for the SPT data release used the lowest 
multipole is I ~ 650. 

In the last panel of Fig. [1] we compare the spectral 
shapes of typical scalar and vector magnetic modes with 
an amplitude of 6 nG and a spectra index of -2.5 with 
the sum of the three astrophysical contributions. We 
note how both astrophysical and PMF contributions are 
comparable in the same range of multipoles. 



IV. RESULTS 

In the present work we derive the constraints on PMF 
performing a combinetd analysis of the WMAP 7 year 
lH, S3 and SPT data following Ref. 

We use the latest WMAP likelihood code 
(version v4pl) and associated data available at 
\http:/ /lambda. gsfc. nasa.gov/\ We modify the WMAP 
likelihood by excluding the temperature bandpowers 
between t = 800 and 1200. We use the SPT data release 
relative to the observation of 790 square degrees of the 
sky at 150 GHz during 2008 and 2009. The data spans 
the £ range from 650 to 3000. In order to decrease the 
correlations between the two data sets we excluded the 
SPT bandpowers for I < 800 and we used WMAP 7 
years data in temperature up to £ — 800. 
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FIG. 1. Comparison of magnetized modes for a field of 
BiMpc = 6 nGauss and ns = —2.5 with the as trop hysical 
contributions from SPT 150 GHz data as in [53 , l45l |. From 
top to bottom. First panel shows the WMAP (blue dots) and 
SPT (yellow dots) data compared with scalar(dotted line) and 
vector(dot-dashed line) magnetic modes. The second panel is 
simply an high I focus of the first. The third panel shows the 
astrophysical contributions: SZ effect (dashed green line), the 
Poissonian term (dotted blue line) , the clustering term (triple 
dot-dashed yellow line) and the solid red line represents the 
sum of the three. In the last panel the magnetic contribu- 
tions including the uncorrelated sum of the two (dashed line) 
is compared with the total astropohysical contribution (solid 
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We develop an extension of CosmoMC [48[ in order to 
compute the Bayesian probability distribution of cosmo- 
logical parameters, including the magnetic ones. In order 
to use the small scale SPT data we introduced the con- 
tribution of astrophysical contaminations following the 
scheme given by the SPT collaboration [32j. We modi- 
fied the code following the indications given in (45j . 

We vary the baryon density w& = fij h 2 , the 
cold dark matter density u> c = il c h 2 (with h being 
iJo/100kms~ 1 Mpc _1 ), the reionisation optical depth r 
(not to be confused with the conformal time t), the ratio 
of the sound horizon to the angular diameter distance at 
decoupling 9, ln(10 10 As), ns and the magnetic parame- 
ters i?iMpc (in units of lOnG) and tib- As priors we use 
[0,10] for BiM P c/(10nG) and [-2.9,3] for n B (> -3 in 
order to avoid infrared divergencies in the PMF EMT cor- 
relators). Together with cosmological and magnetic pa- 
rameters we varied also the parameters describing the as- 
trophysical residual contributions : -Df ^ , -D3000' ^3000 ■ 
We used the prior [0 , 100] for the three astrophysical pa- 
rameters. 

We assume a flat universe, a CMB temperature 
Tomb = 2.725 K and we set the primordial Helium frac- 
tion to j/hc = 0.24. We restrict our analysis to three 
massless neutrinos (a non-vanishing neutrino mass leads 
to a large scale enhancement in the power spectrum of 
CMB anisotropies in the presence of PMF [15| and would 
not change our results) . The pivot scale of the primordial 
scalar was set to fc* = 0.05 Mpc -1 . In order to match 
the data we lensed the primary CMB angular power spec- 
trum using the lensing tool included in CosmoMC, we have 
not considered the lensing of magnetized angular power 
spectrum. We sample the posterior using the Metropolis- 
Hastings algorithm [49] generating four parallel chains 
and imposing a conservative Gelman-Rubin convergence 
criterion [50] of R - 1< 0.01. 

The results of the analysis performed with the combi- 
nation of WMAP 7 and SPT data are shown in Fig. H 
The constraints on cosmological parameters are in agree- 
ment with the one obtained in 32] since as shown in 
(lOj PMF contributions do not modify the constraints on 
standard parameters. Concerning the PMF parameters 
we obtained -Bim P c < 3.5 nG and ns < —0.25 at 95% 
confidence level. 

Contrary to what might be expected from a mere an- 
gular power spectrum point of view, the magnetic param- 
eters are not strongly degenerate with the astrophysical 
ones. In fact in Fig. [3] we show the distribution of the 
astrophysical parameters versus the magnetic ones and 
note that there is no degeneracy between the astrophys- 
ical models and the magnetic modes. This holds for the 
multipole range of the data considered here. 

We note the improvement given by SPT with respect 
to our previous analysis with WMAP 7 and a combi- 
nation of small angular scale data 1 101 which included 
ACBAR 11, BICEP jU and QUaD [3c|. We considered 
ACBAR [2!| data up to £ = 2000. As far as astrophysical 
residual contamination in WMAP 7 and ACBAR are con- 




FIG. 2. Constraints with WMAP 7+SPT data, in the anal- 
ysis we included the six ACDM cosmological parameters the 
two magnetic ones, Bm P c and ns, but also the astrophysical 
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FIG. 3. Constraints with WMAP 7+SPT data for the two 
magnetic parameters, B\ and ns, versus the astrophysical 
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cerned, in Ref. [10[ we relied on the treatment provided 
with the likelihoods which internally account for point 
source residuals, whereas for the TSZ we adopted the 
standard nuisance parameter which represents the am- 
plitude of the TSZ signal, assuming the Komatsu-Seljak 
template based on [5l| scaled accordingly to the WMAP 
and ACBAR frequencies. With that analysis we derived 
the following constraints on the amplitude and spectral 
index of PMF: B 1Mpc < 5.0 nG and n B < -0.12 at 95% 
confidence level. Similar CMB constraints - of the order 
of 6 nG at 95% confidence level - with similar data sets 
were obtained in 
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FIG. 4. Constraints with WMAP 7+SPT data, without the 
inclusion of astrophysical residual contributions. In the anal- 
ysis we included only the six ACDM cosmological parameters 
the two magnetic ones, B\ and ub- Dashed lines are for 
comparison the results of the analysis which includes the as- 
trophysical contribution. 



A. Importance of astrophysical residuals for 
magnetic parameters 

To investigate the importance of the astrophysical con- 
tamination of small scale data for the PMF constraints 
we performed an analysis with the same combination of 
WMAP 7 and SPT data as the previous one but with- 
out taking into account the astrophysical residual con- 
tributions to the angular power spectrum, which means 
setting all the three astrophysical parameters to zero. In 
Fig. [4] we show the results of the analysis, we note how 
thought there is no degeneracy between magnetic and as- 
trophysical parameters the absence of the astrophysical 
contributions in the angular power spectrum results in 
a bias for -BiMpc and tib, which would lead to a tenta- 
tive detection of a few nG amplitude with an ns ~ —1 
spectrum. 



B. Constraints on causal fields 

The results of the analysis with WMAP 7 and SPT 
data shows that positive spectral indices n B > are al- 
lowed only with a very small PMF amplitude. 

In the following we constrain the amplitude for ns = 0, 
2, 3. The results of the analysis are: i?iM P c < 5.6 x 1CT 1 
nGauss for nn = 0, Biu D c < 6.6 x 10~ 3 nGauss for 



tib = 2 and i?iM P c < 7 x 10~ 4 nGauss for ub = 3. 
As expected from the results of the analysis for vari- 
able spectral index the constraints on -BiMpc are much 
stronger than the generic ones. These tight limits will be 
important for their implications on PMF generation and 
evolution. 



V. IMPLICATIONS FOR PLANCK 

In our previous work [Tol j we have analyzed the capa- 
bility of the Planck satellite [52| to constrain the ampli- 
tude of PMF. In Ref. we have performed a MCMC 
analysis with Planck simulated data with vanishing PMF 
obtaining -Bim P c < 2.7 nG at 95% confidence level, a con- 
straint stronger than the present one obtained with the 
combination of WMAP 7 and SPT data. However the 
two constraints are not directly comparable, in fact in 
our analysis of Planck forecasts we did not include any 
contamination by astrophysical contributions of the sim- 
ulated data as it should instead be expected in real data 
on small angular scales. In particular, the presence of ex- 
tragalactic contributions on small angular scales has an 
impact on Planck data and might degrade the constraint 
on PMF amplitude. We now investigate this issue follow- 
ing the treatment of astrophysical contamination, which 
has been developed for Planck data, given in Ref. [36| . In 
Ref. [36| are considered the three main contributions for 
Planck frequencies from 70 to 353 GHz: the Poissonian 
term both for radio and infrared galaxies, the clustering 
term for infrared galaxies and the TSZ. The contributions 
to the angular power spectrum are modelled with empir- 
ical parametrizations based on Planck early data results 
for the point sources and on theoretical models for the 
TSZ. We performed a MCMC analysis with Planck sim- 
ulated data with the combination of five frequencies, 70, 
100, 143, 217, 353 GHz with updated beams and noise 
characteristics from [53[, 54 1. In Fig. [5] we show the 



comparison between the results on PMF amplitude with 
and without astrophysical contamination. The result on 
PMF amplitude without astrophysical contamination is 
-SiMpc < 2.4 nG and is represented by the dashed line 
in Fig. [SJ This result, without astrophysical contamina- 
tion, is updated with respect to our previous work since 
the simulated data are now based on an extended mission 
duration and the sensitivities and angular resolutions are 
updated to the last in flight real performances. The re- 
sults of the case where the astrophysical contamination is 
considered arc represented by the solid line in Fig. [31 we 
note how as expected the constraints on PMF amplitude 
are degraded and in particular are -BiMpc < 3.6 nG at 
the 95% confidence level. 



VI. CONCLUSIONS 

We have studied the constraints on a stochastic back- 
ground of PMF by a combination of WMAP 7 and SPT 
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FIG. 5. Constraints with Planck simulated data with (black 
solid line) and without (red-dashed line) extragalactic resid- 
ual contribution. We followed the treatment for the residual 
contribution as in [36j] for both the simulated data and the 
MCMC analysis. 



data. On the basis of previous works, we have included 
in the analysis only the dominant magnetic contributions 
by scalar and vector magnetized perturbations. In order 
to not introduce biases in the magnetic parameter con- 
straints we have considered the contamination by astro- 
physical residuals of the SPT data. The dominant con- 
tributions are given by unresolved point sources and in 
particular radio and infrared galaxies and by galaxy clus- 
ters. We have considered both Poissonian and clustering 
terms for point sources and the SZ effect for the galaxy 
cluster contribution. To model the contributions to the 
angular power spectrum of the three signals we have used 
the templates provided by the SPT collaboration [lH, [45[ . 
These templates require three additional parameters to 
describe the amplitude of the three astrophysical signals. 
Therefore together with standard six cosmological pa- 
rameters we have included in our analysis not only the 
PMF amplitude and spectral index but also the three 
new parameters for the astrophysical residuals. We per- 
formed a MCMC analysis with the eleven cosmological, 
magnetic and astrophysical parameters and we constrain 
SiMpc < 3.5 nG with preferred positive spectral index at 
95% confidence level. The results do not show any strong 
degeneracy between magnetic and astrophysical parame- 
ters which is compatible with the multipole range of SPT 



data (£ m ax ~ 3000) used. Comparing these results with 
the previous constraints with data by WMAP7, ACBAR, 
QUAD and BICEP [H El, which were of the order of 
-BiMpc < 5 nG, we note a very good improvement for the 
constraint on BiMpc with the use of SPT data. 

We updated the forecasts for Planck that we obtained 
in our previous work considering an extended mis- 
sion duration and the updated in flight performances we 
obtained J5im P c < 2.4 nG. This results do not consider 
the astrophysical contamination of the data which in- 
stead will likeli affects real Planck data. Therefore to 
obtain more realistic forecasts we performed an analysis 
which takes into account extragalactic residual contribu- 
tions to Planck simulated data following the treatment 
in !36| . The results we obtained show a significant degra- 
dation of the constraints on PMF due to the presence of 
extragalactic contributions: -BiMpc < 3.6 nG. The con- 
straint from Planck realistic simulated data is therefore 
compatible with the one we obtained with the combina- 
tion of WMAP and SPT. 

We investigated the impact of astrophysical residuals 
on the magnetic parameters with a MCMC analysis using 
the same data but neglecting all the three astrophysical 
signals. The results show biases in all parameters, both 
the cosmological but most of all in the magnetic ones. In 
particular the PMF amplitude shows a tentative detec- 
tion of few nG. This result obviously is not related with 
the nature of PMF but it is completely due to the astro- 
physical contamination and demonstrates the necessity 
to properly consider possible residuals to constrain PMF 
with small scale CMB data. 

The results presented here confirmed a previously 
noted trend which prefer negative ns- Since ns > is 
mainly related to causal generation mechanism, we have 
shown again how causal fields are allowed with an ampli- 
tude much smaller than the nGauss. 
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